


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1994 


A simple analytical model for asynchronous 
dense WDM/OOK systems 


Huang, Yun-Yao 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/28534 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


+ 
a) 
cai 


a 
i yegt tell 


Ta acee ie 

ule 

pce tee 
erg, 


Ce 
sit dere 4! Gas 

Veen FLere 
ty Les ee 
Lae) (> & 


3 eee 
' Se Uy a ’ 
wilt, “1 “4 3 
> 


6 a 

pe Le 

The tL A we 
etsy, 


a : 
“eng teSyuce F 
” eM 


eg ie "4 
= 
Te | rr aoe ww Se} 


4% 
ee ar) 
rt es 


Pe pee 


< 
-* 


o 3%" 
oteeres 45 
“wes 
tye ‘ ba a er 
fa 
atian a \ dati 


tt Wee pares 


H 
aft Ve, , 
wets Vag, 
x. 


ry < 
WN EE te, 


2. 


—. 34 
esti were, ovee 
ater? 
“« 


: Yon 
at tte, 


an 
: : £otlgie 
%, es . ! \ 


a ef 
Meee eed, 


‘4 
TI epee 
€.4el 


a; 


o 


= 8% @f 


a he mee) 
1+ , 
— eh ay x6 
: . 
mcf ee 24 Tang “t 


« seg 
as = 
rn teby gh 


= 
ie 
t ‘Fg.0 
A ota "hit 
. . mh 25% 
the a pal ge iter et egand, 


ay h 
eet yel'y 
Was 

t 


4 ¢ 
Ul alesse 
. a ee 
takai eas 
a Oe BE t 


it ate C%AY prraae 
‘ . % a a ‘ < ; , Rs 4 ers 
" 4 ete 2 
1 ' : 4 


. . 
oer Hed "le %e te 


«jhe 
2 
tor e% , «oi : LU “re ir a 
dtaonw ’ ¢ ' ; 
nina 1 ; ae PRT LE 
: ¥ Pas] eee tes 
rtrel nse Ags 
ba lis ry 


ote my ys 
Cede gt ae 
oh 


4 ay 
s 
ihe et 


4" 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 








Approved for public release; distribution 1s unlimited. 


A SIMPLE ANALYTICAL MODEL 
FOR ASYNCHRONOUS DENSE WDM/OOK SYSTEMS 


by 
Yun-Yao Huang 


LTCOL, Taiwan Air Force 
B.S, Chung-Cheng Institute of Technology 1978 


Submitted in partial fulfillment 
of the requirements for the degree of 
MASTER OF SCIENCE IN PHYSICS 
from the 


NAVAL POSTGRADUATE SCHOOL 
June, 1994 





REPORT DOCUMENTATION PAGE From Approved OMB No. 0704 


Public reporting burden for this collection of information is estumated to average | hour per response, including the time for reviewing instruction, searching 
existing data sources, gathering and maintaining the data needed , and completing and reviewing the collection of information. Send comments regarding this 
burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, and 
Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503 | 


1. AGENCY USE ONLY 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED | 
(Leave Blank) June 1994 Master's Thesis ) 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 
A SIMPLE ANALYTICAL MODEL 
FOR ASYNCHRONOUS DENSE WDM/OOK SYSTEMS 


6. AUTHOR(S) YUN-YAO, HUANG 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 
Naval Postgraduate School REPORT NUMBER 
Monterey, CA 93943-5000 


9. SPONSORING/MONITORING AGENCY NAME(S) AND 10. SPONSORING/MONITORING 
ADDRESS(ES) AGENCY REPORT NUMBER 










11. SUPPLEMENTARY NOTES 
The views expressed in this thesis are those of the author and do not reflect the official policy or position of the 
Department of Defense or the U.S. Goverment. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. *A 


13. ABSTRACT (maximum 200 words) 

We derive the closed form expression for the bit error probability of asynchronous dense WDM systems employing an 
external OOK modulator. Our model is based upon a close approximation of the optical Fabry-Perot filter in the receiver as 
a single-pole RC filter for the signals that are bandlimited to a frequency band approximately equal to one sixtieth of the 
Fabry-Perot filter's free spectral range. Our model can handle bit rates up to 2.5 Gb/s for a free spectral range of 3800 GHz 
and up to 5 Gb/s when the power penalty is | dB or less. 


14. SUBJECT TERMS 15. NUMBER OF PAGES 
. 79 
16. PRICE CODE 


17. SECURITY CLASSIF-|18. SECURITY CLASSIF- |19. SECURITY CLASSIFIC- |20. LIMITATION OF 
ICATION REPORT ICATION OF THIS ATION OF ABSTRACT ABSTRACT 








PAGE 
Unclassified Unclassified Unclassified UL 


NSN 7540-0 1-280-5500 Standard Form 298 (Rev. 2-89) 





Prescribed by ANSI Std. 239-18 


ABSTRACT 


We derive the closed form expression for the bit error probability of asynchronous 
dense WDM systems employing an external OOK modulator. Our model is based upon a 
close approximation of the optical Fabry-Perot filter in the receiver as a single-pole RC 
filter for signals that are bandlimited to a frequency band approximately equal to one 
sixtieth of the Fabry-Perot filter's free spectral range. Our model can handle bit rates up to 
2.5 Gb/s for a free spectral range of 3800 GHz and up to 5 Gb/s when the power penalty 


is 1 dB or less. 
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I. INTRODUCTION 


Asynchronous wavelength division multiplexing (WDM) systems have been 
increasingly proposed as an attractive alternative to coherent optical frequency division 
multiplexing (FDM) systems [Ref. 1-5]. Although asynchronous WDM systems with 
direct detection do not have the channel capacity of coherent optical FDM systems, they 
are much less costly to implement. Furthermore, present filter technology enables the 
designers to tightly pack the channel, resulting in asynchronous dense WDM systems that 
can provide aggregate bit rates of many terabits per second (1 Tb/s = 10” b/s). 
Asynchronous dense WDM systems are particularly attractive in the area of undersea 
surveillance where hundreds of sensors and data collection sites are envisioned being 
merged onto single-fiber superhighways through massive data fusion. Other applications 
call for relatively low-bandwidth data collection over many months to be dumped quickly 
to a remote recording site in a matter of minutes. This "collection-and-dumped" 
compression can demand total data rates on the order of hundreds of Gb/s. Long distance 
between data collection sites and a remote recording site requires the use of optical 
amplifiers. Therefore, it is necessary to pack all channels within the optical amplifier 
bandwidth. 

An asynchronous dense WDM receiver with on-off-keying (OOK) modulation can 


be modeled as shown in Fig. 1. Conceptually, the analysis involves two main operations: 


1) a convolution operation to evaluate the signal at the output of the optical filter, a 
Fabry-Perot (FP) filter in our investigation, and 2) the integration of the output of the 
photodetector. Evaluation of bit error probability by the numerical analysis of these two 
operations has been carried out in [Ref. 6], with a number of approximations made to 
reduce the computational complexity. In this investigation the FP filter is shown to be 
well approximated by an RC filter within the frequency range | 7-7, |< #SR/ 20x, where 
FSR is the free spectral range of the FP filter [Ref. 7] and f, is the FP filter center 
frequency. For example, given /SR = 3800 GHz, the approximation works very well for 
| f -f, | < 60.5 GHz; that is, the effects of adjacent channels within a 121 GHz bandwidth 
centered at f, must be included, while all others can be neglected. This simple model 
agrees well with [Ref. 6] as demonstrated in Section III. Furthermore, this model enables 
us to obtain a closed form analytical expression for the bit error probability for which 
numerical results can be obtained with little effort. Our investigation shows that this 
simple model provides accurate results as compared to those in [Ref. 6] for bit rates up to 
2.5 Gb/s when the effects of two adjacent channels are included with FSR = 3800 GHz. 
Actually, when the power penalty relative to single channel operation is 1 dB or less, there 
is virtually no difference in the effect of four or two adjacent channels. Thus, for this 
power penalty criterion, this simple model can handle bit rates up to 5 Gb/s for a FP filter's 
FSR = 3800 GHz. 

In Section II the closed form expression for the decision variable, and consequently, 


the bit error probability assuming all channels are bit asynchronous as in [Ref. 6] is 


derived. Section III presents the numerical results which include the bit error probability 
versus the signal-to-noise ratio as function of the FP filter bandwidth and channel spacing, 
and the power penalty (relative to single channel operation without filtering or with 
filtering but no intersymbol interference) versus the channel spacing as a function of the 


bandwidth. Finally, a summary of results appears in Section IV. 


ll. ANALYSIS 


The receiver model for the asynchronous dense WDM system is shown in Fig. 1. 
The desired signal is filtered by a Fabry-Perot (FP) filter that rejects adjacent channels. 
The photodetector is assumed to have a responsivity & (A/W). The detected current is 
amplified by a low noise amplifier that contributes a postdetection thermal noise n(/) with 
spectral density N, (A*/Hz). The decision variable at the output of the integration is 
compared to a threshold & to determine whether a bit zero or bit one was present. 
A. INPUT SIGNAL 

For convenience, we designate channel 0 as the desired channel, and channel & as 
an adjacent channel where & = -M/2, ..., -1, 1, ..., W/2 with M an even integer. We 


consider the equivalent lowpass (complex envelope) data signal in channel 0 and channel k 


as follows: 
0 
bo(t) = 2 bo; Pr(t— iT) (1) 
: Wx ({—Tx) (2) 
b,(t) = 2 bx. je! IU?) Po(t — (IT + Tz)) 
where 


7: bit duration 


5, ,€ {0,1}: bit in channel 0 in the time interval (77, (/+1)7) 


b,, € {0,e/%} is the /” bit in channel & in the time interval (/7 + tz, (7+ 1)T+ tx) 


, : a phase offset between channel & and channel 0 and is assumed to be uniformly 
distributed in (0, 27) radians 


@, : radian frequency spacing between channel & and channel 0 with o, = -@, 
Tt, : 1S the time delay between channel & and channel 0 and is assumed to be 


uniformly distributed in (0,T). 


The function P,{¢ - i7) is defined as 


el Sire eT 
rane ie) — Lo. otherwise (3) 
In both (1) and (2), the non-negative integers Z, and LZ represent the number of bits in 
channel 0 and &, respectively, that proceed the detected bits b,,. The received 


asynchronous dense WDM equivalent lowpass signal at the input of the FP filter is given 


by 


M2 
r() = VP $0 VP b(t) @) 





where P 1s the received optical power. 











B. FABRY-PEROT FILTERED OUTPUT SIGNAL 
The FP filter can be characterized by the following equivalent lowpass transfer 
function [Ref.1,7] 


l—p I—A—p 
A(T) = rani ° =p 


I-p Ap 





7 ay... . 2m, iz (5) 
l-poos(ae}ipsin(zs) —'-P 

where p is the power reflectivity, A is the power absorption loss (zero for an ideal FP 
filter) and FSR 1s the free spectral range. For | | < FSR/20 and assuming A = 0, we can 


approximate H (/) as follow: 


Hf) = aa = ee 


(=p Hye 1+ Go) FSR 
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= aw , (fl <FSR/20n (6a) 
where 
FSR(1- 
C= ---— (6b) 


The free spectral range FSR can be related to the full width at half maximum (FWHM) 


bandwidth B and the finesse F’ of the FP filter as 


te 
el 


HSI 





oa (7) 








Thus if the signal is bandlimited to | f| < FSR/20n, we can truly approximate (5) with a 


single-pole RC filter with the following transfer function and impulse response 





A f)=—s (8) 
l+7—— 
“iy=ce~ , t>0 (9) 


Figures 2a-b show the magnitude and phase (radians) of H( 7) of the FP filter in (5) and 
its single-pole RC filter approximation given in (8) forp = 0.99, F = 312.6, B= 12.16 
GHz and FSR = 3800 GHz. Note that as the frequency increases, the phases of the FP 
filter and the RC filter differ markedly, but the magnitudes of their transfer functions 
remain identical and attenuate rapidly. When | 7 | > FSR/20n, the magnitude of H(/) is 
very small, and therefore, the effect of adjacent channel interference beyond this frequency 
range is negligible. Fig. 3 shows the normalized impulse response of both FP and 
single-pole RC filters. In summary, the above approximation is valid for asynchronous 
dense WDM analysis when the filter finess F is large or equivalently the FWHM 
bandwidth Bis small since the equivalent lowpass signal must be bandlimited to about 
ey | = FSR/20n. 

This approximation has been used in [Ref. 5] to study spectral efficiency of optical 


FDM/ASK systems, which involves the evaluation of the decision variable for worst-case 


analysis using the eye diagram technique. Since we are interested in the detected bit 5 ,, 


in the time interval ( 0,7 ), we consider the output filtered signal s(#),0<t<T given by 


s(t) =sa()+sis(O+sacit), O<tsT (10) 


where 
S, (t): desired signal 
5), (¢ ): intersymbol interference (ISI) signal 
S 4c) (t): adjacent channel interference (ACI) signal 
These signals are evaluated using (4) and (9) as follows (detailed derivations in Appendix 


A): 


sp(t) = VP bo,0 h(t—A)dar 


= J/Pboo(l-e“), O<t<T (11) 


(i+1)T 


= 
si(Q)}=VP X boi J ht-A)da 


= : 
=/Pe* Y bo ifeT el) QersT (12) 


i=—-Lo 


Saci(t) 


M2 Za (/+1)7+t; . 
= JP elie De; | h(t- Aj/elouA—te) qr ree th 
k= Mi 2 l=-L IT+t; 


f 
+by-y | A(t-Aje*OWa_ 


O<t<t 
—T+t, 
tz . 
+ by-y J A(t — rjeloe Wd %<t<T 
= k 
f 
= Jo (A-Tx) 
+ bx AC Ae aM | ry tial 


M/2 —2 1 pt ‘ ; 
=/P > [2 x0, Tbx, evet & : = e(etjo iT +i) (9 (ctjou)T ==) 
k=-MI2 I=-L lj 





a; e IKK : : in 
+n(0, th) bE-} ee ee (e(cvex = @ (C47 x) I+t%)) 
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+n(tx, D)bz-1 & Tmo is: ial eae 
Cc 


+n(TE, T) bro eet ar (e (c+j@x)t _ e (CO x)te)) (13) 


0<t<T 


hise= > |. . 
where 7 (¢1,t2) = 06 otherwise iS 2 unit pulse between 1, and 1, 


The FP filtered output s(t) 1s detected by the photodetector which produces a 
current of @ | s(t) |? Amps. This current plus additive white postdetection thermal noise 
current from the amplifier is integrated by the integrator to obtain a decision variable for 
the threshold detector. 

C. DECISION VARIABLES 
The decision variable Y appearing at the integrator output consists of the signal 


component X and noise component V 


Y=X+N (14) 
where 
if 
X=) als(nlat (15) 
0 
ip 
N =] n(t)dt (16) 


We note that N is a zero mean Gaussian random variable with variance N,/. Substituting 
(10)-(13) into (15) we obtain the signal component X as a function of the three parameters 
tx, C7, and w,7, which represent the effect of intersymbol interference and adjacent 
channel interference. 
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D. BIT ERROR PROBABILITY 


For a detection threshold & and an ISI/ACI bit pattern b = {b,,, 5), $; 1 =-Ly, ...-1 
= -L, ..., 0; k= -M/2, .... M/2, k#0; the conditional bit error probability of the OOK 
signal represented by the Gaussian random variable Y in (14)-(16) is given by [Ref. 8] 


a-—X0 (5) 


P.(b) = Oe) + 20) (18) 


where Q (X) is defined as 


& 








U(x) = 





(19) 


os 


and X, and X, are the values of X in (17) for 5,, = 0 and 5,, = 1, respectively. The 
average bit error probability P_ is obtained by taking the expected value of P, (db) in (18) 
over all bit patterns 6. The minimum bit error probability is obtained by optimizing over 


the threshold a. In summary, given a, we calculate the following expectations: 
ee =E {Pe(b)} (20) 


Pe, min =minE {P.(b)} 
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HI. NUMERICAL RESULTS 


In this section we present the numerical results for a) bit error probability versus 
signal-to-noise ratio Z = RPJ7/No as a function of the normalized channel spacing 


(normalized to the bit rate) 
— OF _ 
f= oe MAS) (22) 


where Af= @,/ 27k is the equal channel spacing inHz, andb) power penalty versus 
normalized channel spacing / as a function of the FP filter parameter cT=(nx/,/p )BT 


where B 1s the filter full width at half maximum bandwidth(FWHM). For our model to be 
valid we set M = 2, that is, we constraint the signal to be bandlimited to within two 
adjacent channels, thus we incorporate the degradation wv the two nearest adjacent 
channels. We observe that for bit rates of 1 Gb/s or less, ou: model is valid up to M = 10 


and very little difference is observed between If=2 and M=10. Also, we observe that 
there is little difference between 14=2 and M=4 when [210 for bit rates up to 3 Gb/s. 


In all results we set L, =2 and L=1. 
Figures 4-5 show the minimum bit error probability versus signal-to-noise ratio 


( Z ) for c7 = 5 and 10, respectively. Along with each curve, we also show that of the 
synchronous case, i,e, tT =0,L, = 2 and L = 0, and that of a single channel (SC) 


operation without filtering or with filtering but without /S7. In Fig. 4 we observe that a 


large degradation occurs due to JS/ for cT = 5 which represents a narrowband filter. As 


the FP filter bandwidth is made larger as in Fig. 5 with cT = 10, the /S/ is reduced but the 
ACT increases. 
In our model, we are constrained to M = 2 for the case under consideration. We 


use an FP filter with FWHM B=12.16 GHz, free spectral range FSR = 3800 GHz, finess 
F=n /p /(1—p) = FSR/B = 312.6, and c = 38.4 GHz; 1/T = 2.56 Gb/s, and cT = 15. If 


the power penalty related to single channel operation 1s to be less than | dB, then by 
Fig 6, a minimum normalized channel spacing of / = 12 must be used. Equivalently 
from [Eq. (22)], we could use Af = //T = 12*2.56 GHz = 30.4 GHz, (i.e. the channel 
spacing is a multple of the bit rate). In this model the farthest adajcent channel for M= 2 
is twice the channel spacing which is 60.8 GHz. This verifies the assumption | f-f, | < 
FSR/20n = 60.5 GHz, where f, is the FP filter center frequency. This result agrees well 
with that in [Ref. 6; Figs. 6,9, M/F =0.4, a = 0,2]. Thus we incorporate the degradation 
caused by the two nearest adjacent channels. We observe that for bit rates of | Gb/s or 
less, our model is valid up to M = 10, and very little difference 1s observed between M = 


2 and M = 10. Also, we observe that there is little difference between M = 2 and M =4 
when / 2 10 for bit rates up to 3 Gb/s. In all results we set L,=2 andL=1. 
Figure 6 also shows the power penalty for an asynchronous dense WDM system 


relative to a single channel operation at the minimum bit error probability of 10°’. This is 
the required additional signal power (dBW) for the asynchronous dense WDM system to 
be able to operate at the 10°” bit error probability achieved in the single channel system 
with a SNR=12dB. The asynchronous dense WDM system is /S/-limited at 2.15 dB, 0.95 


dB, 0.6 dB, and 0.5 dB in power penalty for cT = 5, 10, 15, and 20, respectively. It is 


eee ee — 








seen that for a 2.3 dB power penalty, the normalized channel spacing can be as close as / 


= 6 (i.e., a channel spacing of six times the bit rate) for cT. = 5. If the power penalty 
criterion is 1 dB, the normalized channel spacing 1s / = 12 for cT = 10, 15, 20. We remark 


that although the exact transfer functions of the FP filter is used in [Ref. 6], a number of 
approximation have been made to obtain numerical results. The approximations are 1) 
the /S/ is obtained by modeling FP filter as a single-pole RC filter [Ref. 6, Eqs. (4) and 
(36)], 2) approximating the finite integration with an infinite integration in the calculation 
of ACT [Ref. 6, Eq. (15)], and 3) the beat interference is ignored. On the other hand, the 
[ST and ACT in our investigation are obtained by modeling the FP filter as a single-pole RC 
filter, using finite integration and including the beat interference. Since the results in our 
investigation and in [Ref. 6] agree well, we conclude that approximations are quite valid. 
We also note that our results also agree well with the simulation carried out in [Ref. 1, 
mae 17). 

The above numerical results shown in Figs. 4-6 are obtained with an optimized 


threshold setting. Figure 7 shows the power penalty for fixed threshold a = @P7/2 which 


is the same optimum threshold for single channel operation (midpoint between the 
received power for bit zero and bit one). It is seen that the performance of an 
asynchronous dense WDM system is quite sensitive to « for a narrow band filter. An 
additional 1.5 dB is observed for c7 = 5 for J > 8, and 0.4 dB for cT = 10 for J > 12. 


Negligible degradation is observed for cT = 15, 20 for J > 16. 
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Figures 8-9 show the power penalty versus normalized channel spacing as a function 
of FP filter parameter c7 for the worst-case analysis with optimal threshold and fixed 
threshold, respectively. The worst-case bit pattern is fixed to produce the minmum X, 
and maximum X, where X, and X, are the values of X in Appendix-A equation (4-21) 
with 5, , = 1 and 4,, = 0 respectively. 

We observe that the power penalty for the worst-case analysis 1s only slightly larger 
than that of the exact analysis for / > 10 shown in Fig. 6. Similarly the power penalty for 
the worst-case analysis with fixed threshold is only slightly larger than that of the exact 
analysis with fixed threshold for J > 10 shown 1s Fig. 7. The reason for this is that for 
large channel spacing (J > 10), the ACT effect is small, so the AC/ bit pattern has a small 
influence on the power penalty. 


Figure 10 shows the normalized optimal threshold for the exact analysis shown in 
Fig. 6. It is observed that a ~0.4 for J > 10. Note that the normalized optimal threshold 


for the single channel operation is a = 0.5. 
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IV. CONCLUSIONS 


We have presented a simple model for the analysis of asynchronous dense WDM 
systems employing an external OOK modulator. The only approximation that we use 


involves the modeling of the Fabry-Perot filter by a single-pole RC filter assuming the 
equivalent lowpass signal is bandlimited to the frequency range |f| < FSR/20x. This 


model enables us to obtain a closed form expression for the bit error probability which 
previously could only be obtained via numerical analysis [Ref. 6]. For an FP filter with an 
FSR around 3800 GHz, our model can include the ACI effects of two adjacent channels 
for bit rates up to 2.5 Gb/s. Our numerical results show that this model agrees well with 


that in [Ref. 6]. 
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Figure 2: Spectral characteristics of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 
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Figure 3: Normalized rimpulse response of the Fabry-Perot filter and the approximated 
single-pole RC lowpass filter 
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Figure 4: Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for c7=5 
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Figure 5: Probability of bit error versus signal-to-noise ratio as a function of normalized 
channel spacing for c7=10 
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Figure 6: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter cT 





Power Penalty (dB) 
+ on 


G) 





Figure 7: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter c7 with a fixed threshold a=0.5 
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Figure 8: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter c7 for worst-case with optimal threshold a=(X,+.X,)/2 
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Figure 9: Power penalty versus normalized channel spacing as a function of Fabry-Perot 
filter parameter c7 for worst-case with fixed threshold a=0.5 
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Figure 10: The normalized optimal threshold versus normalized channel spacing as a 
function of Fabry-Perot filter parameter c7 
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APPENDIX _A 


Derivation of Formula 





IL INTRODUCTION 


Desired channel =: Channel 0. 


Adjacent channel : Channel k. 
k =-M/2,...-1, 1,...MW/2, M: even. 


Bit in channel 0 in 7, time interval (/7, (i+ 1)T) : 6,,€ (9, 1} 

Detected bit in channel 0 in (0,T): 6,,€ {0,1} 

Bit in channel k in/, time interval (/7'+ ty, (7+ 1)7+T,) : 5, pelo tw) 
by € {0,2}, j= J-1 
@, frequency spacing between channel & and channel 0, © = —O-x. 


;, ; phase offset between channel é and channel 0, uniformly 
distributed between (0, 27). 


T, : time delay between channel & and channel 0, 0<t;,<T. 


0 
Data signal in channel 0: bo(f) = >, bo; Pr(t—iT) (et) 


i=—Lo 


0 
Data signal in channel k: ,(#) = >) by jel) Pr(t— (IT + ty) (1-2) 


l=-L 
L,, L : positive integers. 


_ 1,0<t<7T 
je T (t) a { 0 otherwise 


The received signal at the input of the FP filter of channel 0 1s : 
M/2 
r(t)=JPbo(t}+ >» VP b(t) 
k=—M/2 
k#0 


P: received optical power. 
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The output of the FP filter 1s: 


ro(t) = i W(t —d)r)dr 


= 


- JP Pf h(t—boad.+/P Sf Ae A)bAAAA 


k— Mie —0o 


= VP boo J h(t— A)Pr()dh. 


+ JP & boy [ h@—A)Pr(—iTyan 


I=—L 0 —OO 


WP 2 _{ b> fe J h(t del Prd — (IT + ty) dd 


+by-} i h(t —A)el#O-W P(A. — (-T + t,)) dr 


—00 
oro) 


thio | h(t Ae Pr(A-—te)dr F (1- 


—oo 


h(t) is the equivalent lowpass impulse response of the FP filter of channel 0. 
Since the detection interval is 0<t<T7, we only need to evaluate 


S(t) =r, (), O<t<T andwrite S(t) as: 
SID = S30 + Sis + Sac Y Daisey 
S; () : Desired signal in channel 0, the Ist term (/=0) of (1-3) 
Sis; () : Inter-symbol interference, the 2nd term (i < -/) of (1-3) 
Sic{Y: Adjacent channel interference, the 3rd ierm (>) of (1-3) 


All evaluated in 0 <t< T. 
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I. CALCULATE S, (t) and S,5/(t) 


In this section we are to calculate the easy part of S(t) --- S, () and S,., (). 


A. FP filter: 


The lowpass equivalent transfer function of FP filter is 








HG l-p 1-—A-p l-p 1-—A-p 
7 janfFSR* |J-p 2 a 
(me I-p 1-pcos(=4)+jp sin(=) I-p 
p : power reflectivity. 
A: power absorption loss ( 0 for ideal filter). 


FSR : free spectral range. 


Since {<< Fsr (for operating frequency range), we can approximate 
H (f) as (assume A =(0): 








_ | ee 
H@) = ] me Lif = 2nfp ] 27f 
(—-p)P ise = lao esr ce 
FSR(1- 
where C= ae 


for FP filter, we also have : 


FSR _ =r 
B 





, B: full width at half maximum bandwidth (FWHM) 
of half power bandwidth. 


Using inverse Fourier transform, we find 


ee eee, 
hy = oe PCAN se 
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B. 


Calculate Si (t): 


S,() =AP boo | h(t-2) Pr(a) dr 


—c&O 


f 
= JP boo | h(t—2) da 
0 


f 
= JP boo | ce a 
0 
f 


= JP boo ce a4 dn 
0 


= JP boo(I- e~et) Q<f=R (2- 


(a) For Bit"l" : S,@ = V/P(l-e“) 0<t<T 
(b) For Bit) Oss eae 0<t<T 


Calculate_S,., (t) 


5S. =vP 3 bos J Me- h) PRCT 


i=-Lo 
G+1)T 


Sj 

=JP X bo | h(t-A)aa 
i=—Lo iT 
_] (+1)T 


=/P > bo; | SND Gi 


i=—Lo iT 


= 
Pe dD bo; tM! — e®*) dr awh (2-1 
i=—Lo 


a5 


Ill. CALCULATE S(t) 


In this section we are to find S,., with the procedure : 


A. When /< -2 

B. When /=- 

C. When /=0 

D. Summary for S,., () 
A. Whenl<-2 





Integrate the first term of 5, (t) in (1-3), we have : 


Py ¥ bi J h(t —2) eO-W P(A, — (IT + 14) dr 
ei 2 l=-L 


9 (/+1)7+t; 
=JPYY by | AG-A)elO-wWan, 
kl IT+t, 
5 (+1) T+te 
=/P>> bx) | cel) ej hte) I, 
ko 


IT+t, 


-) (+1) T+t 
=JP YY bi | ceretiann) eletionkgy 


ko! IT+t, 


=2 
-_ JP » » eee =(Ci+j@;T;) e (ctjo «)(IT+7,) (eee kT _ 1) 
k 1 


= pp »y Sy bx eet .— @ (470TH) (9 (CAO WT _ 1) O<1<T (3-1) 
(i ltjz 


c 


=! 
=P > ie. e~ AR) (glerjoNlTFT) _ o(ctjoxIT ) O<?t<7 = (3-1a) 
k 


ie Gel 
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B. 


when l= -! 
Bel [=-1, 0<t<T 


Integrate the second term of 0, (t) in (1-3): 


M/2 22 
JPY bya | Alt-A) WP 7(X.— (-T+ tA 
k=-—M/2 
k#0 


—0o 
f 


= JP Dobe | h(t-A) ean 
k 


—T+t, 


f 
= Ae » bp-} | Ce (CHJ@ kT) e(CHOx)A 
k 


—T+t, 


=e » bp_1 en"! eo (een — eletjox)(-T +t) 
k ; ly 


O<t<T ( 
= JP » eel Cs = e Ait) e(c7ox)F) 
y lye O<t<T , (: 
Bolg — |. ea fa] 
M2 & | 
JPY bi-a { A(t-A) ef Pr — (-T + ty))da 
k=-M/2 Be 
k#0 
Tk 
= JPY bp | h(t-2) eWay 
k —7T+t, 
= —ct @TPKk (ctiox)te(1 _ p-(chHon)T 
UP 2 bi-1e Tee (1-e ) 1 <i<T 
ips , | 
= aye » aH ect) de e~(ctiax)T) Tt. < fae (3 
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when l=0 
C-1. 1=0, 0<t<T, 


Oo 


VP Y bio | Ae—Aje!™* Pr(A-Ty) = 0 O0<t<t (3-4) 
k me 


wee, 1 —0,7, <r 


Integrate the third term of 5, () in (1-3): 


oo 


JP 2, bro | A(t- V)eFox(A—ti) Pr(A- Tan 
k —co 


t 
= JP » bro | h(t- A)eIou(Ate) dh 
k Tk 








JOTy 
= Pp b eet € ~ eg (Ct7Ox)i —e (CHO )Tx 
* ” oh : mt — (3-5) 
EJPY bk (e@x(-%) pelt) Tr<t<T (3-5a) 
, lis ; 
Summa or_s iol 
Mee: U(t,,t,) —U,-U, i, t= 0, ti} Sto 


where U. (t) = me a c20 


is the unit step function. 
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then 7 (t, t,) = { 


l, f1StSto 
0, otherwise 


is a unit pulse between ¢, and f,. 


hence for 0<t<7 ~~ S,.,() can be written as: 





M/2 
Se (t) = JP SS > 70(0, T)bx: pt) aa e ae e COMETH) (p (CHOK)T aa 1) (3 
ly 


(3-6) 


al 2 /=-L 


M? 
+J/P > n(0,t,)bg-1 eo —S ee: a (e (cA7K) 9 (Cty@x)(-T#%%)) 
k=—MI2 +] 


M/2 
JOpTy 


[VP XS mt, Ne em 


i=l DE 


e (CH Mi TE ( ] —eevowl ) 


M/2 


—j® : P 
se > (Tr, T)b x0 oh ae tk = (e (C+jWx)t — e(eteute)) 
k=—M/2 le 


=7, +2, 42.6 Zi oe 


ZZ. Ge) 
Tig (3-2) 


T 
Zy-1—> (3.3) 


Zs (3-5) 


k,0 


ge, 


6 
B- 


(3. 


a ee er ee ———_ 


a 


IV. WORST-CASE ANALYSIS 


Worst-case for S(t). 

Worst-case for & | | S(t) |? dt - output of integrator. 
Special worst-case analysis. 

Special worst-case at the output of the integrator. 


SA Bp 


he 


Worst-case for S(t) 


A-l| Worst-case, IS]: 6,,=b., L,=% 


From (2-2), we have : 


-] 
Se) = JP b_ eet! — 1) > eicl 








= JP be“(eF—1) D (e“?)! | 
h=l h=-i 
= JP b_e~'(e°T — 1) er 
| gr) | 
= JP b_e~"(e? —1) 
= JP b_e~* O<t<7 (4-1) 
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A-2 Worst-case ACI: (2p =i, i =e 


M2 | - | 
Sw.(t)=JPb YX {n0, NE (eVT — Ip Y (evo yy 
=—MI2 ae = 





] ee ea ail 
+J/P b py 1(0, Te) oe (eJ Out T) _ p(T) p (c47i)T) 1 





+JP b > 5g Gop. T) ] (eC) _ ect) eran)! ) 
k 


ieee 


+ 1(Te, Na (eF x(t) _ g-ll—T)) 


= JP b > (0, Nx Bes) pe Cine ' 
k Te : 


a [p 5 > (0 .:)—— (eJ Oct) a elitr) e (+70)! ) 
2 Or 
k l+j— 





= (2/0) — ecli-tk) e(CHjMg)L ) 
loge 


+VP b> (tz, 7) 
k 


= JPb »y 70(0, Nm (el) eter!) 
k 7e7 


+.)/ P b »2 (0 tN (J Ox) —_ ecl(iHtk) e{e7@x)1) | 
3 Op 
k lie 


= VPbY nO, 0<t<T 
k Ta 


4] 
| 


4G 
B. Worst-case for & IS(@O|? at - output of integrator 
0 


Output of photodetector : Z |S (t) |’ 


ISO = (SO + Sift) + Suc, (0 | 


= |S,OV + SiO + |Sacr Y | + 2REES, (Y Sisy(0 } 


ip 2Re{ Sz (t) Such 3 ee {Sis (0 Sach j 


e 
breed 


From (2-1) 


Sp’ () F =|S, (|? = Pb, ,? (1-2e% + e"), 0<1<T 


‘a 
| ISa()l2dt = PTB} of - 21 - e~*) +L - e*7)] 
B-2 From (4-1) 


a) | = Pb? e**, 0<1t<T. 


T 

we b2 —2C 
[ Is@lar == (1 —e72°7), 
0 
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(4-3) 


(4-4) 


(4-5) 


(4-6) 


(4-7) 


B-3 From (4-2) 


M2 MB al 
See AON = Pb >a 7. 
pen Viper Viv, (a Uli): 

k#0 m#0 


ell@p-@m)t 6 —}(@ pt p-O mtm) 
= P\b|? >») > —————— 0<1<T 
rom (Ute) (7) 


i 
fis "“cy(O)l2dt = Plal2 » eo e —}(@ ,t,-@ mtm) J eker nda 
q om (ly) (I) 


eH e Kone Onin) Gea vy == 
= x = e m —. l 5 k m 
_ - 2 HES y 4) Ig) ( ) 
PT\bl2 2 





Tebe 2 = 
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2Re {Sz(t) Sig (D} = 2Phoob_(e~! — ee") Q<reT (4-11) 
f PTbo0b 

[ 2Re {Sa(t) SiG (O} dt = Par —eP)? von 
0 


B-5 

2Re {S*(0) S%*.,(1)} = 2PhooRe {b (1 - ey } 

= 2PhooRe {b > a (eon — eXCFOM)\ Qe t<T a) 
k ghar 


2PTb @ pT pT —(c—j@ ,.)T : 
—** Re (bo ( ae eek’-1 _ l-e ——) (4-14) 


ey 2 





T 
| 2Re {S3() Sue(O}dt = 
ww 


2Re {Sig () Sie} = erik 5 ~(c-jou)t 
(Sisr OD Sac} = 2Pb_Re (od ye Pere? 15) 





2Re {S¥* (1) SM. (t)\dt === Re (b 2 a: (l eee eS iis 
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9 
a | 


C-2 


i 
ho 


C-3 


C-5 


Special worst-case analysis 


O, = 27kl/T, I: positive integer 
Tk = O 

bo; aS {0,1} 

bx} = bro =be 10, el?) 

Re{bse {0, cos } 

lbI2 € {0,1} 


T 
[ \Se(pl2dt = PT b) (1-20 -e*7) +E -e *)] 
0 


T 
silos b2 = 
[ Is¥s@l2.dt= = a —e7) 
0 
in (4-9) 
if 
| ol(@r-Om)t ip = pee: 





ISacr(t)| Seat = 


= od- 


ip 
| 2Re {Sa(OSig (0 } scdt = —2= (1 -e*7)? 
0 


Gi 


T 

Re i? 2PTbyoRe{b}  _, 
[ 2Re {Sp(DS4ey() Jace = POO (eT - NE a 
0 





vi 
Ges 2PTb Re{b ee 
#] 2Re {Sisi (OSie(0 }scdt = ASN —e“D a 
0 
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(4-5) 


} 


(4-17 


(4-18) 


(4-12) 


(4-19) 


(4-20) 


D. Special worst-case at the output of the integrator 


T 
yur = & [lsvolzar 
0 


= ¥ {all terms in section C]. 


HPT § bo p1 — ey ey 


2a 
_ +25 (1- er") 
— (l-e<")? 
+ (lol? +22 Re {O}(b-—boolID Taare} (4-21) 
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V. EXACT ANALYSIS - GENERAL CASE 


Now, we are to encounter the most complicated situation, for clarity, we 
first list the construction talbe, for which we are to follow - 








S(t) 
a 
| | | 
Sicil t) Sisil( t) Sal t) 
| | | 
Zi, Zo 2 Table 5-1 
_ 
| | 
Z%,0 ae —| 
IS@l? 
| 
F E _ D_ Aae B A 


| | | | | | 
2Re{Sis(OSsc(O}  2WRe{SR(HS4c)(t)} 2Re{Sa()Sis(O} — |Suer(0)|? ISis()l2 Sa? 
| 


| 
C-6 C-5 C-4 C-3_C-2_ C1 
| 


| | | 
2Re{Zo'Zi,} 2Re{ZoZ0} rRe¢zizity Ize l2 Ize? [zz 2 
| 
| 


omre 3-2 C35 
| | | 
Tt 
le 5- 2Re{Z7_2Zu5} IZ2,,\- Zoe 


We first calculate each term in Table (5-2), then integrate it. For each term 
belongs to ISucr(O)|?, we also discuss two cases: W,=W,, and Wr#Wm . 
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Since the computation is very complicated, here we just write down 
the results for each term, and omit the intermediate steps. 


ii 
A. {|Ss(ol2dt= PT; [1 -A + Se (4-5) 
0 


B. From (2-2) 


= 
Sis)? = Pe2"[ bo (et DT — eT)? ear I (5-1) 
i=—Lo 
T =| 
| ISisr(O|2 dt = rang = erly 2 bo s(etVT — eieTy)? (5-2) 
0 I=-L0 


C. Calcualte |S fll - refer to (3-6) 


Cl 


=: ae 
IZp I? = Pe'x(0, TI DY by Ee ec HOW(T+W) (g(evO0T _ 1)|2 (5.3) 
kl) We 


T 
yy . 
a8 9) bal PT ‘ JO;,t, : : 
| | Z| Ay — ST 4 2cT)| py » bx ‘att e CAO) ITH) (9 (CH @4)T . 1)| 2 (5-4) 
0 c 
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C-2 T, =n eee 


or (ge (CHI@4)! _ (CHM: N-T414)) 
-Oon ) 


1% 
(+j)1 


IZ) (= P SOS nOwer sae 
k om 


(e (C—j@m)t =e (CJOm\— ET mM 


—}( m) ; 
CHOI OMEM Fie Om)t — ACH On) g (CHO N-T4) 





=P >> 0, eo ae 
kom 


(l4j—Ey 12) 
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i) K=M,O, =On, 


a =o eis [c p+ = (e7 (C17 ute ” Lye OWE te) 
lH)? Nae 


+ (e —(cj@K)tk — | Ve (c—jo ;,)(—T+tx) 
oF hares 


—s(e = 1) e (C+J@ «)(—T+tx) e (CJOrwy-L +Tx) 
va lbx-11? 1 2ct pyeie 
oF cy OEE OI 
1 Reta (een? — plevorN-71))\) (5-6b) 


0-3 3 terms for |Z" | 


C-3-1 Ty = max(Tx, Tm) 


e HOt, -@mtm) 


vas 2_—p —Jch 7 - ———— 
| ti | a 2 Mtg, i») € bx. 1Dm-\ el) 


[elevoute(] — eXcvonl) ele Jom)tm(] — ee Jom) (5-7) 


fl 22Par= fz Plat: 


(a)k#AmM, Or FOm 


Dre 1b; = =5ye5 ICT ge J@Rt~—-@ mtm) 
Bee 6-2) 
_—* Yue SIZ) | 


[e (c+jo te] aa e (ct wl) e(cJo m)tm (1 -_ ec jo my) (5-8a) 
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(D) k=, Op =o, 





PT lb, 11? = = E ; 
= aaa rome — er( ele) | +e 2cT _ Jo CF Re fel!) 


C-3-2 


e Io Kte-Omtm) 
ace) etic 


VA aes © a 2d, W(tg, T) bi,0bm07 
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TKO 
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l+j-3 
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ise ar 


l mS = A ee 
=a 2cl _p *ate) e (CTO x)TK po (C JO m)Em) 


Ok =On= or, 


et Jol? 1 2(-cT+ct;) 
= of 24 eb [cel -—cte+ ail —e k)) 
+2Re{ ewes (e~ CHOW _ p-(erjox)tey\] 


te 


5] 


(5-8b) _ 





(5-9) 


(5-10a) 


(5-10b) 





iz gt tp —-W mtn) 


DRe{Z: Ee = 2PRe{ 2, a Wty T)b,- \bmon, Fy] om) 


[en(cvo mt p(x) (| — g-(et7o0)T) 


224 elecion)in elerionte(] — gv) 5.41) 


ar 
[2Re{Z2_ ONG J eZ Zixo}at , fo kezAmandk=m 
0 


PRetD x by, 1b ae : e Charli tei gle eee!) e (CHO K)Te 
(Sea Ss) 


= e —(C+J@ m)tg _ e —(C+/O »)T ) 





[ 


7: 


—(e-**s _ ya) e (JO m)tm 7} (5-12) 


a2 


C-4 


—}(® -T,.-@ mtm) 
2Re{ ZZ 1 =2PREAYD > n(0. tm)b,1b%,_. ———— 
ko om l=-L aire 


[ e AG Onn) e (c+jo x)U1+tx) (e (C470 )T _ ] ) 


_, cr e (c-j@ »,)(-T+t m) e (c+jo ;)(UT+t x) ( e (c+7@x)T _ | )] } 5. 


i Tm * 
| 2Re{Z§Z¢ Jat = | 2Re{ZpZo" de for k#m and k=m 
0 0 


Ba e HO xtp-@ mtm) . ‘ 
~ Re » b SO e (ct7@ fl | e (cto n)(T+t,) 
ar ap 19m ae & ) 





_p—(ct+/om)tm _p72etm) ; 
(eS ae ge (C-J0 m\(-T+tm) J (5-1: 
1+j-4 : 
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C-5 


* (@,-T.-@ mtm) . | 
DRe{Z) VAS = 2PRetd 2 5 (Tm, Lp es (ection) _ Le Cen)lT+14) 
l ica) ice 


[bx1b,, a | e2clp(cjo m)tm (| _ e(CJOm)!) 


+bz 1b), 0 e(C+J@ m)t 


—bijbji9 eee om ]} (5-15) 


fe T 
[2Re{ZjZh ydt= | 2Re{ZiZ" at 
0 


Tm 


UE R fd 5; * as all rod (cio _ Jp (e47@ NIT +x) 
bom tet dees 


(e726 aa Bagel (A = e (ce Danae m)Tm 


[ 


yale ey 
2 





+ eae ie —(c+j@ m)Tm ae e (ca m)t) 
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C-6 


DREW aan The product is only survived when Tm <f< Tx 


e HO ¢t,-o@mtm) 


= DERC T (Tm, Tk - = 
2 a Come te) 8m 
[Dp-1 be, ene TO uM (CF mitm (| — eC Om) 
by 1b* _,e-2etelHON-T#4H) 9 (CF m)tm(] — gC m)T) 
Oe bin. el (@ k-Om)t 


= k-1 bio e(cJo i) (ej m)Tm 





Si] b 0 e (C4J® m)t (C47 x)(—14 tk) 
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(5-17) 
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‘g Tk 
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of REALL Te ee —. 
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lj 
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l+j— 
reksla (e72etm = e~ 26% (CAFO MT 47K) @(CJO m)tm J} } (5-18) 
§ ig 
D. | 2Re{Sp(t) Sis} dt = [ 2Sp(0 Sisr(t) at 
0 0 
= = 2Pbho0 5 bo; eet 0H) fe eu) dt 
i=-Lo 
= ma bo o(1 = = 5" by; (Orne = gaa (5-19) 


i=—Lo 
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E.  2Re{Sz() Sac} 


etPRR 
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= 2PRe{ZE m(0, T)(e~ — ebook. 





x (C (c+jo,)T _ | Je (c+jo ;)(/T+t,) 
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+2 bo,0by,- re: Gora ek Nie (a) -——_ ae 





c 
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0 


= aye y bo j (eT — eT) Re TS: by(1-e7 a 


Ye (ctjo,)T _ 1) e(cryo KT p Gur 
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ae a Ee [; — [eo - e JOR e(CJonl) — x(e-# — ectke2cl) J } 
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VIL EXACT ANALYSIS - SPECIAL CASE 





Now, we put two constraints into all terms of exact analysis to form a 
special case ; 


l.@,==, I: positive integer. 
2.7% =T for all adjacent channel k. 


For the purpose of normalization for future use, we set Q = 


S14 


so that O.T, = OT = aT = 2nkla 


| 
Still, we follow Table 5-2. : 





ates ~T) (l-e7267) 




















cl+j2nkl 


i =] | 
| Sis()Pat = FE. —e 7) [ Le doe"? — el) 7? 

0 i=-Lo (6-2) 
C. Special case for ACI 
| 
Cl 
fof PT es “pala : 
| IZaileae _ = ee eter & = by" male" = jl} elcT o acl j2nklo 2 | 

- =, ~ly= 

0 
_ eter (e2aer — Beale lg 1)? | 3 bye? | 2 (6-3) 


3 


C-2 


fizslear= fizilear 
0 0 


(a)kK AMO, z Om 


deat OF, 
— 2 ee 
=(cl)*P Td 2 (cT+j2nkl)(cT- “faa 
]—e V2ma(k-m) 
| J2nI(k-m) 
ecl_p(a-IecT+j2nmla 


cT+j2aml 


eel _p(a-I)cT—j2nklo. 
cT—j2nkl 


e2(a-IcT_p-2cT 
gee -eut (6-4a) 


(b) k = m, Ox = Om 


lo, I? oe) 
¥ y) k-1 SS 
=({c/) PY ds prone (cT)?-{(2nkl)* [a+ Zoe 
—cT_,(a-1)cT+j2nkla 
ut ecco gee — 
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C-3 3 terms for |Z! |? 
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ie 3 
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(6-8a) 


(6-8b) 


(6-9) 


(6-10) 


if 
[ 2Re{ Sa Sis) hat 
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el -D | | 
= a (i - Cay 5 bo, (Ceo =e) (6-11) 
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=2PI Re \ a : 2(cT+j2Tkl) 
3 boobr-1 cT (1-e a) (e pee taste Eee eee Val 
cT+j2nkl | j2nkl cT-j2nk] 2cT 
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Re 2 2 2(cT+j2nkl) 
5 ye. icT ie (Caeser) (cen enone!) 
= cT- yoael | elle  ¢T-j2nkIl 7 2cT 
> b k-1 Pa saa, (l-e~*?) 
2(c7T+j2tkI) 
kee by ocT (Os (eet _g(e-2c7) (6- 13) 
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APPENDIX _B 


Programs 





% exact analysis with optimal threshold signal out of the integrator, the signal X contains desired 
signal and ACI & ISI and postdetect noise ,The formula for the BER 1s : 


——- _ p(b) 


M(L+1)4L 
DAS oak aay pattern 


where p(b) =; (>- )@(=)™ 


ae i a—X'o(b_14/2---O.072 T-M2---TW2) 
[ ie | Co 75 - Abdur at_yn...atun 
MO MO 


J Nor 


2m ik 
X1 (by... LMI? «++ — 
= [...J] ie,,| ONS db van..dbuan At_yn...atmn ] 
MO MO 0 


M is the number of adjacent channels. Here we assume all channels have the same phase @ and 


time delay T i.e, Ox = 0, Tk = T_ hence for our model the power of (1/270) and (1/T) are fixed 
to 1, therefore we only have two integrals and two arguments. 


M=4; k=[-M/2:-1 1:M/2]; 

m=k; 

% produce the controlled matrix b to control 64 different bit patterns 
m1=[ zeros(1,32) ones(1,32) ]; 

m2=[ zeros(1,16) ones(1,16) zeros(1,16) ones(1,16)]; 


m3=[]; 
for 1=1:4 
m3=[m3 [zeros(1,8) ones(1,8)]]; 
end 
m4=[]; 
for 1=1:8 
m4=(m4 [zeros(1,4) ones(1,4)]]; 
end 
m5=[]; 
for 1=1:16 
m5=[m5 [zeros(1,2) ones(1,2)]]; 
end 
m6=[]; 
for 1=1:32 
m6=[m6 zeros(1,1) ones(1,1)]; 
end 
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b=[m1;m2;m3;m4;m5;m6]'; 


% signal to noise ratio range in dB 
RPTN_DB=[10:0.5:25]; %RPTN_dB 
ppp=10.*(0.1*RPTN_ DB); 
lenl=length(ppp); 


% solalph function is provided by Professor Randy L. Borchardt 


n=10; 
[bpx,wfx]=grule(n); Yobpx=bpy ,wfx=wfy 


% single channel 
BERO=0.5*erfc(ppp/8%0.5); 


pp=[]; thresh=[]; thresh is not normalized threshold 


for CT=[5 10 15 20] %cT is Fabry-Perot filter parameter 
if CT==5 
qqq=[2:12 17 20]; 
elseif CT==10 
qqq=[3:12 17 20]; 
elseif CT==15 
qqq=[4:12 17 20}; 


elseif CT==20 
qqq=[5:12 17 20]; 
end 


pe=[]; thresh]=[]; 


for I=qqq 
BER=[]; thresh2=[]; 


for RPTN=ppp 
ap=linspace(0,182,11); | % approximated thresholds 


x3=0; % first few loops to find out the threshold which make the 
% BER minimun don't care about the scale 


for 1=1:64 x3=x3+... 


solalph(‘xx00',0,2*pi,2,bpx, w£x,0, 1,2,bpx, wfx,ap,b,CT,1,..k,m,RPTN)-+-... 
solalph(‘xx1 1',0,2* pi,2,bpx,wfx,0, 1,2,bpx, wfx,ap,b,CT,1,.k,m,RPTN); 


end 
[val,ind]=min(x3); 
lef=ap(ind)- 16.6; 
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if lef<O ,lef=0; end % to aviod the threshold go beyond the negative side 


ap=linspace(lef,ap(ind)+16.6,11); 
% two more time to find alpha 
x3=0; 


fori=1:64 x3=x3+... 
solalph('xx00',0,2* pi,2,bpx,wfx,0, 1,2,bpx,wfx,ap,b,CT,1, .k,m,RPTN)+... 
solalph('xx] 1',0,2*pi,2,bpx,wfx,0, 1,2,bpx, wfx,ap,b,CT,1,.k,m,RPTN); 
end 


[val,ind}=min(x3); 
lef=ap(ind)-3.1; 


if lef<O ,lef=0; end 


ap=linspace(lef,ap(ind)+3.1,11); 
% three more time to find alpha 
x3=0; 


for i=1:64 
x3=x3+... 
solalph(‘xx00',0,2*pi,2,bpx,wfx,0, 1,2,bpx,wfx,ap,b,CT,1, |k,m,RPTN)+... 
solalph('xx1 1',0,2*pi,2,bpx,wfx,0, 1,2,bpx,wfx,ap,b,CT,1,..k,m,RPTN); 
end 


[ val,ind}=min(x3); 
lef=ap(ind)-0.57, 


if lef<O ,lef=0; end 


ap=linspace(lef,ap(ind)+0.57,8); 
% four more time to find alpha 
x3=0; 


for 1=1:64 
x3=x3+... 
solalph('xx00',0,2*pi,2,bpx, wfx,0, 1,2,bpx,wfx,ap,b,CT,i,.k,m,RPTN)--... 
solalph(‘xx1 1',0,2*pi,2,bpx, wfx,0, 1,2,bpx,wfx,ap,b,CT,i,|k,m,RPTN); 
end 


[ val,ind ]=min(x3); 
ap=ap(ind), 
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% after find optimal alpha use double integration 
qq3=0; 


for 1=1:64 
qq3=qq3+( dbgquadm('xx00',0,2*pi,2,bpx, wfx,0, 1,2... 
bpx, wfx,ap,b,CT,1,.k,m,RPTN)+... 
dbgquadm('xx1 1',0,2*pi,2,bpx,wfx,0, 1,2.... 
bpx,wfx,ap,b,CT,1,..k,m,RPTN) )/(256*pi); 
end 


BER=[BER qq3]; thresh2=[thresh2 ap]; 


if qq3<10%(-15) %for save time only interesting in ]0“(-15) 
li=find(ppp==RPTN); 
BER(ii+1 :len] )=5*10’(-116)*ones( 1,length(ii+1:len1) ); 
thresh2(iit+ 1 :len])=(ap+2)*ones(1 ,length(i+1 :len1) ); 
break 
end 
end 


pe=[pe;BER]; threshl=[thresh] ;thresh2]; . 
end 


pp=[pp;pe]; thresh=[thresh;thresh] ]; 
end 


time2=toc; 
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r 


a. o-_ aan aaa 


= —— 


% worse case form appendix setting optimal threshold equal (x0+x1)/2 
M=4: 
I=linspace(0,6, 101); % set I value in linspceinteger to find minimum I 
ee N DB=0:0.2:20; % x-axis dB range 
RPTN=10.%(0.1*RPTN_ DB); % change to ratio 
% single channel 
BERO=0.5*erfc(RPTN/8“%0.5); 
%find optimal alpha 
pp=[]: 
for CT=[5 10 15 20] 
x3=0; 
for k=1:0.5*M 
pe —X3t2./(1+(2* pi*k* I/CT).“%2); 
end 
x0=(1-exp(-2*CT))/(2*CT)+x3*(1+2*(1-exp(-CT))/CT); 
x1=1]-2*(1-exp(-CT))/CT+(] -exp(-2*CT))/(2*CT); 
[val,ind ]=min(abs(x0-x1)); 
J=ceil(I(ind)); 
J=(I(ind)); % our minimum I value 


% set different I and to find the BER 
ss=I+]:1+10; 
BER=[]; 
for J=ss 
x3=0; 
for k=1:0.5*M 
X3=x3+2/(1+(2* pi*k*I/CT)“2); 
end 
x0=(1-exp(-2*CT))/(2* CT)+x3 *(1+2*(1-exp(-CT))/CT); 
x1=1-2*(1-exp(-CT))/CT+(1-exp(-2*CT))/(2*CT); 
BER=([BER;0.5*erfc( RPTN*(x1-x0)/8%0.5 ) J; 
end 


pp=[pp;BER]; 
end 


semilogy(RPTN_DB,BERO.'--'RPTN_DB,pp(:,:)) 


axis({10 19 10%(-15) 1]) 
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